Helical edge resistance introduced by charge puddles 
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We study the influence of electron puddles created by doping of a 2D topological insulator on its 
helical edge conductance. A single puddle is modeled by a quantum dot tunnel-coupled to the helical 
edge. It may lead to significant inelastic backscattering within the edge because of the long electron 
dwelling time in the dot. We find the resulting correction to the perfect edge conductance. Gener- 
ahzing to multiple puddles, we assess the dependence of the helical edge resistance on temperature 
and doping level, and compare it with recent experimental data. 



The realization that a boundary separating a 
topologically-nontrivial insulator from a conventional one 
should carry delocalized electron states P, Q has led to 
the prediction of such states in concrete materials and 
their experimental observation ■ One of the stunning 
theoretical predictions is that in 2D the zero-temperature 
electron transport along an edge is reflectionless, as long 
as time-reversal symmetry is not broken, which should 
lead to the quantization of the edge conductance [2:| . 

Experiments with HgTe quantum wells of the ap- 
propriate thickness confirmed the existence of highly- 
conducting channels in a nominally insulating state of 
a heterostructure [^-Q. The Fermi energy Ep in a het- 
erostructure was tuned by a gate to reside in the gap 
between the valence and conduction bands. The values 
of conductance G measured under these conditions were 
indeed close to the predicted quantized value Go = e^/h 
per edge, but only for small ~ 1 x 1 /im^ samples. Devi- 
ations AG = Go — G towards lower conductance values 
were clearly seen in larger samples [6l-[T]|. 

In short samples, AG fluctuated with gate voltage. 
The temperature dependence of G has not been system- 
atically measured yet, but the existing data indicate it 
to be rather weak. These observations should be con- 
trasted with the theoretical predictions of a strong tem- 
perature dependence of electron inelastic backscattering 
rate, with a characteristic scale set by the band gap Eg- 
Depending on the model, AG scales as oc [T/Eaf' or 
(T/Ea)'^, unless time-reversal symmetry is broken [3.IT2I- 
Il4l |. Spontaneous symmetry breaking is improbable for 
weak electron-electron interaction (noting the high di- 
electric constant, k « 13 i5|) and in the absence of a 
high density of magnetic impurities. Luttinger liquid ef- 
fects d, [13, 0, UHl are also suppressed in short samples. 

The existing theory considers inelastic electron 
backscattering by either uniform interactions alon g an 
edge 0, m, Ell, or at isolated points 0, El, El, 113 ■ 
Helical edges formed in a semiconductor heterostructure 
are likely to deviate considerably from either limit. The 
structures are doped [HI, El, EH ; the presence of charged 
donors and acceptors results in a non-uniform potential 
landscape for electrons. These inhomogeneities are not 
point-like because of the long-range of the Coulomb po- 
tential. Moreover, the topologically non-trivial insula- 
tors are in fact narrow-gap semiconductors with a typical 
gap of only Eq — lOmeV for HgTe quantum wells fi- 



ll. El, EH, Ell- To place Ep inside the band gap, an ap- 
propriate gate voltage is applied, so that the gate charge 
balances out the uncompensated donor (n^) or acceptor 
(ua) charge density. The joint effect of the gate and ion- 
ized dopant atoms may lead to the formation of electron 
and hole puddles in the quantum well, cf. Fig. [U similar 
to the known phenomenon in compensated bulk semi- 
conductors (19i] . In the existing measurements doping 
varied from Ud ^ 3.5 • 10^^ cm^^ to ria ^ 5 • 10^" cm~^ 
El, EH , and the results indeed seem to indicate that 
a lower doping level improves the quality of the edge 
conductance quantization. Furthermore, the uncovered 
strong sensitivity of the edge conductance to the poten- 
tial of a scanning probe [10[ may imply the presence of 
puddles, i.e., spontaneously formed quantum dots, in the 
vicinity of the edge. 

In this Letter we elucidate the role of tunneling be- 
tween an edge and a quantum dot on the edge conduc- 
tance. Elastic processes involving electron dwelling in the 
dot do not lead to any backscattering. However, dwelling 
enhances the inelastic backscattering by increasing the 
time electrons interact with each other. At tempera- 
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FIG. 1. (Color online) Electrons moving along a helical edge 
tunnel in and out of puddles created by the inhomogeneous 
charge distribution in the heterostructure. In the puddles 
electrons may undergo inelastic backscattering. The main 
contribution comes from puddles whose distance d from the 
edge is within a strip where the resulting level width F ~ 
T, cf. Eqs. (O-®. The strip width is the tunneling length 
A = v/Eg- Summation over the puddles yields the average 
resistivity, Eq. (fT2|l . 
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tures T < 6, the dwelling-time effect makes the conduc- 
tance correction strongly dependent on the position of 
the Fermi level Ep with respect to the dot energy levels, 
and on the tunneling widths F of these levels {5 ^ Ea is 
the mean level spacing in the dot). At a given temper- 
ature T, the tallest peaks AG?''^'" oc (T/S)^ in AG{Ef) 
are produced by levels with F ~ T, see Eq. ([5]) and Fig.[T] 
Such peaks in AG{Ep) are of widths ~ T, and the "peak- 
to-valley" ratio is ~ (S/T f. 

Dots, or puddles of charge carriers in a quantum well, 
are formed accidentally by fluctuations in the donor den- 
sity dll, [l^l- We establish a crossover value tiq of Ud 
below which puddles are rare. At <^ uq the density 
of puddles. Up, is exponentially small in no/ud- In short 

— 1/2 

samples of length L < Up only a few puddles are in 
the vicinity of the edge, resulting in mesoscopic fluctua- 
tions of G with the gate voltage. This model agrees with 
the results of scanning-gate experiments [lOj] and could 
explain the variations of G with the back gate voltage in 
earlier experiments (^-Q , if the condition Ud ^ tlq would 

— 1/2 

hold there. For longer samples, L ^ Up , many pud- 
dles couple effectively to the edge. That leads to edge 
resistance, R oc npL{T / S)^ , which varies smoothly with 
the gate voltage and possibly greatly exceeds the quan- 
tized value /i/e^. At the same time, the "bulk" hopping 
conductivity, which is proportional to factors exponen- 

— 1/2 

tially small in Up and T/S, may still remain negligi- 
ble. In this case, current would flow along the edges, 
despite edge resistance being high compared to /i/e^, as 
observed in Ref. [ll|. The model would also explain the 
earlier measurements on larger samples, if the con- 
dition Hd < no would be satisfied [our crude estimate of 
no, Eq. pUj) . turns out to be too low for that]. 

We start by considering a helical edge coupled to a 
single quantum dot via a point contact. In the absence of 
interactions, the corresponding Hamiltonian takes form: 



The inelastic backscattering reduces the steady-state 
current I = Iq — A/ from its ideal value Iq = GqV by 
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+ E ^"4^^7(0) +H-C.. (1) 



Here vp is the helical edge velocity, 7 = ±1 = i?, L la- 
bels the right- and left-movers, respectively, and n labels 
the discrete energy levels in the dot, measured from Ep. 
The dot is coupled to the edge at x = by a set of tun- 
neling amplitudes i„. The Kramers degeneracy of each 
discrete energy level n gave us the freedom to pick the 
corresponding eigenfunctions \n'y) in such a way that the 
left- and right-movers are coupled to two different com- 
ponents of each doublet. There is thus no backscattering 
in the free-electron problem. Interaction in the dot, 

^ ^ 2 ^ y ^ni7in272;n373n474''ni7i'^n272''"474''n373' 

,7i 

may lead to inelastic backscattering (hereinafter, we as- 
sume U respects time- reversal symmetry). 
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{El, E2; E^, E4)S{Ei + E2 — E3 — E4) 

-f^,iE,)f^,{E4)il-UiEi)){l-f\AE2))] . (3) 

Here V is the source-drain voltage, AA'^^j^2;-y3^4 = (71 -|- 
72 — 73 — 74)72 counts the net number of right-movers 
scattered into left-movers, f^y^iE) = l/le^^+^^^^f^y^ + 
1] is the Fermi function shifted by ±6^/2, and 
'^7i72;7374(-^ii E2; E^, E4) is the cross section for the two- 
electron scattering process |-Ei7i, -E272) — ^ \E3j3, E4J4} 
between exact left- and right-propagating eigenstates of 
the Hamiltonian ([1]). In general, 5* allows for backscat- 
tering of one {RR — > LR) or two {RR — LL) electrons. 
There are two respective contributions, AGi and AG2, 
to the conductance G = Go — AGi — AG2. 
In the Born approximation the cross section is 
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Here Qnin2 (E) is the noninteracting retarded Green func- 
tion of an electron in the dot. All interaction matrix ele- 
ments must be small compared to F to allow the pertur- 
bative treatment at arbitrary position of the Fermi level 
with respect to the dot levels. This condition is more 
easily satisfied for the off-diagonal matrix elements [l^ 
entering explicitly in Eq. ([4]), than for the diagonal ones 
C^ni7in272;n.i7in272 ^ Ec- The introduced charging en- 
ergy Ec is small, Ec <C 5, if the spacer between the 
quantum well and gate is thinner than the Debye radius 
for electrons in the well. In the opposite case of Ec ^ 6, 
Coulomb blockade may develop. We first treat the en- 
tire interaction perturbatively, and later point out how 
Coulomb blockade modifies the results. We will also see 
that backscattering is dominated by puddles with F ~ T; 
thus Kondo correlations [l3| setting in at the much lower 
temperature Tk F can be ignored. 

Using properties of the interaction matrix elements in 
Eq. it is straightforward to check that in the low- 
temperature limit AGi oc and AG2 oc T^, in agree- 
ment with Refs. [T3- 14 1. For a generic form of strong 
spin-orbit interaction in the dot, all interaction matrix el- 
ements in Eq. ([4]) are of the same order . In this case, 
AG2/AG1 < 1 if T < (5. The proportionality cocflicient 
of the temperature dependence AGi oc T** is a function 
of the dot parameters; in the case of weak tunneling it 
peaks every time a level crosses the Fermi energy. 

Weak tunneling corresponds to small elastic tunnel 
widths F„ = \tn\'^/{2vp) of the levels, F„ < |£„ - e„±i|. 
Then the leading-order approximations for the diagonal 
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and ofF-diagonal (rii 7^ 112) matrix elements of Q^{E) 
read g^^{E ) = [E - e„ + zr„)-^ and g^^^^{E) = 



'ir„2)] \ respec- 



tively. Using this simplification in Eq. Q we find 
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(5) 



for the peak in AGi corresponding to the level £1 cross- 
ing the Fermi level (ei = 0). The peak height and its 
width, \ei — Ef\ ~ Ti, display mesoscopic fluctuations; 
Eq. ^ is apphcable at T < Ti. The peak value l^Gf"^ 
grows with temperature till T reaches a value T ^ Ti. 
At higher temperatures, some of the incoming electrons 
with energies \E\ < T which contribute to AGi are off 
resonance. This leads to a decreasing T-dependence of 
AGf'''' at T>ri, 
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(6) 



and a peak width \ei — Ef\ ^ T. 

In a weakly-disordered dot the Thouless energy Et = 
gS ^ 5 [g ^ 1 \s the dimensionless conductance within 
the dot). The disorder-averaged matrix elements (/7^) of 
interaction present in Eqs. ([SJ and ^ can be evaluated 
using the standard diagrammatic techniques |20|]. Fur- 
ther simplification is possible for the screened Coulomb 
interaction, which is dominated by its universal zero- 
momentum component, leading to ([/^) ^ 5^ / g^ ■ Us- 
ing this estimate in Eqs. ©-(jl]) and dropping numerical 
factors, we arrive at the interpolation 



(AGf"'^) 
Go 



^%(^i^-T) + y^O{T^T) (7) 



for the typical peak conductance as a function of T at 
small charging energy, Ec <C maxjT, F}. 

The backscattering processes leading to Eqs. ([5]) and 
([6|) involve a sequence of virtual states. Those with en- 
ergy deficit |£„| ^ are represented by the denominators 
in the sums over n ^ 1. One of the virtual states, how- 
ever, has two electrons on level n = 1 and brings a large 
factor l/F^ to Eqs. © and ©. It is replaced by l/E}, 
in the presence of charging energy Ec ^ F. For the 
same reason, the cross section Eq. ^ loses sensitivity to 
the energy Ei of one of the two electrons involved. The 
corresponding integration range in Eq. (|3]) is restricted 
then by T rather than F at any T/F. In the important 
(see below) case Ec ~ 5, the two modifications change 
Eq. by a factor - {T/5f ■ max{l,T/F}, leading to: 
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A^-T) + ^^0{T~T). (8) 
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Backscattering in the "valley" (Fermi level in between two 
subsequent dot levels) regime does not involve any low- 
energy virtual state and is not affected qualitatively by 



Ec ^ S. The corresponding estimate, (AG™"'''*')/Go ~ 
T^r-^/g^S^ is smaller than the peak value Eq. ® by a 
factor - (FV(56) . max{T'^,T^}. 

The main contribution to the backscattering correc- 
tion averaged over the position of the Fermi level comes 
from the peak values, Eq. ([5]), as (AG™"°^)/Go is para- 
metrically smaller. Accounting for the peak widths, 
\ei — Ef\ ~ max{F, T}, we find 
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e{T - F) . (9) 



At higher temperatures the above mechanism gives 
way to thermally-activated backscattering processes. 
Those originate only from the diagonal elements Q^^ {E) 
in Eq. ([4]). Since this regime is probably not relevant 
for the interpretation of existing experiments (see below) 
we only sketch the results, deferring a detailed discus- 
sion [231 . There are two types of activated contributions 
to AG. The first one involves transitions within a pair 
of levels, {713,^14} = {711,712}. The other one involves 
more levels, {713,714} 7^ {711,712}, and gains importance 
at higher temperatures (T ^ 5) due to the larger phase 
space available for transitions. At T < 5 backscatter- 
ing is dominated by the two levels closest to Ep, and 
AG - ((5V5^rr)-exp(-£/T) with £ ~ 5. Comparison 
with Eq. shows that activated backscattering exceeds 
^^^poak^ at T > (5/ln((5/F). The distinction between 
peaks and valleys is lost at these temperatures, although 
AG does experience strong mesoscopic fiuctuations at 
T ~ 5 due to the randomness of the activation energy e. 

Now we turn to the typical experimental case 
H, [ll] of a doped, gate-controlled heterostructure. For 
definiteness, we will address the case of 7i-doped sam- 
ples, assuming donors of average density 71^ are ran- 
domly distributed in a plane situated between the gate 
and quantum well; the distances of the donor plane and 
gate from the quantum well are id and ^g, respectively. 
Random distribution of donors creates random poten- 
tial V{y) for the charge carriers in the well. In the ab- 
sence of carriers, the variance of the potential [2l[ is 
{V) = Vi\n{il/[i2£g - £d)id]} with Vo = ^S^eV^ 
(k is the dielectric constant). At the point of full deple- 
tion (the gate charge density is —end) the probability of 
creation of electron and hole puddles depends on the ra- 
tio EG/i2y/(y^). The relation ^/iV^ ^ Eg/2 defines 
a characteristic donor density. 



8nenn{ey[{2£g-ed)id]}' 



(10) 



The carrier puddles are small and rare if 77^ <C 77o ; in the 
opposite limit (77^ ^ uq), puddles are large and separated 
by thin depletion strips. In the following estimates we 
ignore the logarithmic factor in Eq. ([TU]). 

In the limit 71^ ^ tiq, fiuctuations of the bands edges 
with respect to Ep are large compared to Eg/2. That 
allows us to use the linear approximation, £(k) = vk, 
for the electron spectrum in the well j5|. An electron 
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penetrates over a length A ~ v/Eq into a p-n junc- 
tion between the puddles. The junctions are formed by 
spatial fluctuations of the random potential V{y), and 
the typical width of the depletion region in a junction is 
Eg/\VV{y)\. Tunneling is weak if EG/\VV{r)\ > A. To 
estimate the characteristic value of |Vy(r)|, we use [l^] 
the correlation function (|Vy(ri)||Vy(r2)|) V^/\ri - 
r2 P at |ri — r2 1 ^ A and find that the weak-tunneling con- 
dition is Eg/Vq > 1. It is not satisfied at Ud ^ uq; the 
p-n junctions are penetrable, leading to F > (5 and aver- 
age bulk conductivity Cbuik ^ e^/h. Recent analysis p5j | 
suggests that a transition from the topological insulator 
to conductor state occurs at (Tbuik ~ (1.4 — 2.5)e^/ft-. This 
makes the limit ^ tiq unfavorable for the helical edge 
conductance quantization at any temperature. 

In the opposite limit, Ud <^ Hq, the puddles can be 
made rare by tuning the gate voltage. In addition to 
no there is another characteristic density scale, l/cig ~ 
27rQ;'*no, set by the effective Bohr radius as = 2hv/aEG 
[we used m* = Eg/2v'^ for the effective electron mass and 
a — /{kHv) for the interaction parameter]. We restrict 
further consideration to the case a <C 1 and no ^ » 
1/a^, where electrons form a continuous two-dimensional 
liquid at zero gate bias. Gate-induced depletion breaks 
the liquid into puddles, which become sparse once the 
electron chemical potential is shifted by > Vb below the 
conduction band edge |2l| . For a ^ 1 we may describe a 
puddle by the Thomas- Fermi approximation. The puddle 
size w is found by matching its number of electrons N 
and electrochemical potential N/{2m*uP') to the typical 



impurity charge fluctuation N 



1/2 



and potential 



fluctuation Vb, respectively. This leads to w ^ aB and 



,5 ^ a^EG ; 9 - ind/no)'/'/{2n)'^^a . (11) 
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Here we used estimates 6 ^ Vq/N, g ^ v N [20[, and 
assumed w < £g, leading to Ec ^ /{kw). Puddles are 
located at rare strong fluctuations of the potential, and 
are thus far away from each other. Their density n^ is 
estimated as the ratio of the total carrier density to the 

1/2 

number N ~ CLsnJ of electrons in a puddle. To find the 
former quantity we note that the distance to the gate £g 
serves as screening length for the potential fluctuations; 
hence we can divide the sample into roughly independent 
regions of size £g . A region becomes populated by carriers 
only if the local potential experiences an exponentially 
rare fluctuation exceeding Eg/2. The carrier number is 
such that they compensate for the fluctuation [2l|. This 
leads to Hp - l/(^gaij)(nd/no)i/2e-"°/"<*. 

The hopping conductivity facilitated by the puddles 
is proportional to a product of two small parameters: 
the tunneling probability, exponential in — (A^np)"-^/^; 
and the thermal activation probability, exponential in 
—S/T. The latter one remains small at T ^ 5, even 
when approaching the crossover region rid ^ no. Un- 
der the same conditions, the rate of backscattering into 
the helical edge scales as a relatively low power of T/6, 

1/2 

cf. Eqs. ©-([H]). For samples of length L < Up only 



a few puddles occur in the vicinity of the edge. That 
would make AG sensitive to a local probe potential, con- 
sistent with recent scanning gate measurements [13] , and 
may also provide an explanation for the mesoscopic fluc- 
tuations of AG{Ep) in earlier measurements P"!^, [l0| . 
Eq. ([5]) predicts that the largest peaks in AG{Ep) at a 
given temperature T scale as {T /5Y and are produced 

— 1/2 

by levels with F '-^ T. For L ^ Up contributions of 
many puddles add up incoherently, as scattering off each 
puddle is inelastic. The exponential dependence of the 
level widths F esc on the distance d between a 

puddle and the edge leads to a broad distribution of F. 
Hence, by Eq. (|9]), backscattering will be dominated by 
puddles whose distance from the edge is such that F ~ T. 
Summing over puddles (see Fig. [T]), we find that a long 



edge displays resistance R 
resistivity 



PcdgcL with self-averaging 



Pcdgc 



1 1 , 

'^0 9 



(12) 



While podgG (T/S)^, the quantum well hopping con- 
ductivity CTbuik is exponentially small in S/T. There- 
fore, leakage into the bulk at Ud ^ no is insignificant 
for samples shorter than the exponentially-large " leakage 
length" L* = 1/ (cbuikPodgn) , which may explain recent 
scanning SQUID results [ll|. 

Our findings thus match with observations provided 
that Ud < no and T < (5. To estimate 6 for a HgTe/CdTe 
heterostructure, we use Eq. pT|) with Eg = lOmeV and 
a « 0.32 (found with k « 13 and v — 5.5 • 10''cm/sec [H]). 
We arrive at (5 « ImcV, comfortably above ksT in most 
experiments. For the crossover density, Eq. pO)) . we find 
no ~ 3-10^'^cm~^. The doping levels reported in Ref. [Toj 
and in 0-[§] are, respectively, moderate {ua/no ~ 1) and 
high (n^/no ^ 10) with respect to this value. On the 
other hand, from the total resistance of long samples in 
Ref. m we deduce that Cbuik ^ 0.45Go, consistent with 
an insulating bulk [25j . It may mean that our crude 
estimate of uq is off by a factor of 10. The characteristic 
length Up^^^, separating mesoscopic samples from the 
"self-averaging" ones, provides another check. The pre- 
exponential factor in it, ~ lOOnm, is only ^ 10 times 
shorter than 1/xm-long "mesoscopic" samples in Refs. 
Q. That too may indicate that the samples doping was 
close to the true crossover value no. 

To conclude, disorder in a doped heterostructure may 
lead to appreciable backscattering within a helical edge, 
while hopping conductivity in the quantum well remains 
negligible, which may explain some of the recent obser- 
vations p-[i, [13, Eli • The samples doping level n^ appar- 
ently was close to the crossover value no separating the 
regimes of low and high bulk conductivity, as opposed 
to our crude estimate, Eq. PH]) . This complicates anal- 
ysis of the of edge resistance dependence on n^. On the 
other hand, the robust qualitative features of the resis- 
tance T-dependence, in both the mesoscopic and the self- 
averaging regimes, Eqs. ([S])-® and (fT^ . respectively, 
makes its detailed measurement very desirable. 
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